GGPLOTs:

ggplot(df, aes(x=x, y=y))+
  geom_point(aes(col = parcela_bl), size = 3, alpha = 0.5)+
  scale_color_manual(values = eqm)+
  ggtitle("Fall Trees Location")+
  theme_light()+
  coord_fixed()
## Warning: Removed 5 rows containing missing values (geom_point).

ggplot(df, aes(x=x, y=y))+
  geom_point(aes(col = succ), size = 3, alpha = 0.5)+
  scale_color_manual(values = eqm2)+
  ggtitle("Fall Trees Succecional Condition")+
  theme_light()+
  coord_fixed()
## Warning: Removed 6 rows containing missing values (geom_point).

ggplot(df, aes(x=x, y=y))+
  stat_density2d(geom="tile", aes(fill = ..density..), contour = FALSE) + 
  stat_density_2d(col = "black")+
  scale_fill_viridis_c()+
  ggtitle("Fall Trees Density")+
  theme_light()+
  coord_fixed()
## Warning: Removed 5 rows containing non-finite values (stat_density2d).
## Warning: Removed 5 rows containing non-finite values (stat_density2d).

Interactive Plots:

int = ggplot(df, aes(x=x, y=y))+
  geom_point(aes(col = succ), size = 3, alpha = 0.5)+
  scale_color_manual(values = eqm2)+
  ggtitle("Fall Trees Succecional Condition")+
  theme_light()+
  coord_fixed()

ggplotly(int)

Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.